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The amplitudes for photoproduction of two pseudoscalars on a nucleon are expanded in the overall 
cm. frame in a model independent way with respect to the contribution of the final state partial 
wave of total angular momentum J and its projection on the normal to the plane spanned by the 
momenta of the final particles. The expansion coefficients which are analogues to the multipole 
amplitudes for single meson photoproduction contain the complete information about the reaction 
dynamics. Results of an explicit evaluation are presented for the moments Wj m of the inclusive 
angular distribution of an incident photon beam with respect to the cm. coordinate system defined 
by the final particles taking photoproduction of and iY°n as an example. 
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I. INTRODUCTION 



The study of multiple meson production is essential for understanding the properties of baryonic resonances, es- 
pecially of those having sizeable inelasticities and for which only a weak evidence from elastic ttN scattering exists. 
According to the quark model calculation of [l|, at least below 2 GeV some of these resonances must be strongly 
coupled to ttttN and ttt]N channels. Therefore, present experiments on 7T7t and ttt] photoproduction have become a 
center of attention in programs discussed at various research centers, and a number of new accurate data have already 
been reported 0411 ■ 

Improvements in the quality of the data have made it possible to perform rather detailed theoretical analyses of 
photoproduction of two pseudoscalars. The 7T7t as well as irr] models have already been the object of several studies [Tol — 
fl8j . Mainly they cover the second and third resonance regions describing with varying degrees of success the existing 
data and predicting the results of new measurements. A typical analysis is based on an isobar model approach. Its 
key assumption is that the amplitude is a coherent sum of background and resonances usually parametrized in terms 
of effective Lagrangeans. As a rule, the resonance part contains s-channel resonances decaying into ttttN or nr]N via 
intermediate formation of meson-nucleon and meson-meson isobars. As adjustable parameters one usually takes the 
masses and partial decay widths of the resonances as well as their electromagnetic coupling constants. 

Within this method, angular momentum decomposition of the amplitude is ruled by partial wave transitions of the 
resonance states to quasi-two-body states, like ttA or rjA. Clearly, such an approach can not be viewed as a general 
partial wave analysis, since it crucially depends on the assumptions about the production mechanism. Therefore, 
it results in various uncertainties, primarily in the non-uniqueness of existing solutions, since the same observables 
may equally well be described with different sets of parameters. Consequently, in spite of a general qualitative 
agreement among the models, significant quantitative discrepancies still remain. Additional limitations may arise 
from inadequacies of the isobar model description, such as violation of unitarity, nonrelativistic dynamics etc., whose 
impact on the description of the processes under discussion remains unknown. For instance, unitarity conditions may 
be important in the region where many production channels are open. 

It is worth to note that one of the main reasons for the lack of a rigorous partial wave analysis for 7T7t and irr] 
photoproduction is that there is no general recipe to deal with reactions involving three particles in the final state. In 
contrast to single meson photoproduction one faces here the technical problems associated with three-body kinematics, 
where the particle energies and angles are distributed continuously. As a consequence, a conventional partial wave 
decomposition of the final state does not provide a multipole representation for practical applications, primarily since 
there exists a variety of ways to successively couple angular momenta of the participating particles to a total angular 
momentum. 

In this paper we present an alternative method by using a partial wave expansion for the photon induced produc- 
tion of two pseudoscalars on a nucleon, which should be of minimal model dependence. It is based on the correct 
determination of the partial wave amplitudes for these reactions with no built-in prejudices concerning the production 
mechanism. Similar method have been used to analyse pion production in ttN collisions, see, for example Refs. [l9l.[20|. 

The paper is organized as follows. In the next section we introduce the partial wave expansion and construct 
the transition amplitude for photoproduction of two pseudoscalar mesons. In Sect. IIIII we use the so far developed 
formalism to discuss some gross features of ir°ir° and ir°r] photoproduction. Finally, some general conclusions are 
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drawn in Sect. IIV1 



II. THE FORMALISM 



In this section, we collect the formulas used in the present analysis. As a starting point the formal results of 
Ref. [21] are used. There the formal expressions for the helicity amplitudes as well as for the cross section and the 
recoil polarization were derived, including various polarization asymmetries with respect to polarized photons and 
nucleons. 

A. The T matrix 

We consider here the photoproduction of two pseudoscalar mesons, denoted mi and mi with masses M\ and M2, 
respectively. Firstly we determine the T-matrix elements of the electromagnetic m,\mi production current J 7mi7Tl2 
between the initial nucleon and the final mim2-/V state. The four-momenta of incoming photon, outgoing mesons, 
initial and final nucleons are denoted by (w 7 ,fc), (u>i,q\), (w2,<f2), (Ei,pi), and (E,p), respectively. The helicities 
of photon and initial and final nucleons are denoted by A, [A, and v, respectively. In a general frame the transition 
matrix element is given by 

T v \n = Q> v I c\ • J imim2 (0)\Pi,n) > (1) 

where for the description of the final state we choose the final nucleon momentum p = (p,d p ,<p p ) and the relative 
momentum of the two mesons q = \{qi — 9*2 ) = (<7, S q , (f> q ). For the following formal considerations the knowledge of 

the specific form of the current J 7mi „ l2 is not needed. 

After separation of the overall cm. -motion the general form of the T-matrix is given by 

T v \,j. = q, v \Jym im *,\(k)\n) ■ (2) 

It is convenient to introduce a partial wave decomposition of the outgoing final state according to 

(-) <£f, P, v\ = — %h ( l P Q \ u \ip v ) (j P m p l q m q \JM) Dl» mp (<f> p , -0 p , -<j) p ) 

Ipjp m p l q m q JM 

x 2^,(6,, -0,,-^)<-><gp; {{l p \)j p l q )JM\, (3) 



where the "hat" symbol means, for example, l q — y/2l q + 1. Furthermore, l q and m q denote total angular momentum 
and projection, respectively, of the two mesons, l p , j p , and m p orbital and total nucleon angular momentum and 
its projection, respectively, and J and M the total angular momentum of the partial wave and its projection. All 
projections refer to a quantization axes to be determined later. For the rotation matrices D , m we follow the 
convention of Rose [22j . 

The multipole decomposition of the current reads with k — (fc, 9 1 , <f>j) 



J imim2 Ak)^-V^Ti L LO^ L (k)Di lLX (^,8 7 ,-^), (4) 



LM L 

where Oj^ L contains the transverse electric and magnetic multipoles 

o\L = E k + XM m l • (5) 

For the initial nucleon state we have 

= (-!)!+* \\™>) D ll-M^-^)- ( 6 ) 

m=±l/2 

Using the Wigner-Eckart theorem and the sum rule for rotation matrices 



M L m 



J L ^^D)l 2 (R)D L MhX {R) = {-lf-»- M ( ,/ x f 3 )D{ fX _JR), (7) 
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FIG. 1: Definition of the coordinate system in the cm. system. 



one obtains 



Ll v j v m v l q m q JM 



X i -m) (/a a i M )(^;((^^) J ll° AL H^ 

x ^„ lp (0 p , -0 P ) D^ m< (0„ -0„ -0,) D J MX _^,6^, -0 T ) . (8) 
Parity conservation results in the following symmetry relation 

lU-A-^n,, q 7 ) = (-i)*-"-"^^,, n p , n 7 ) , (9) 

where for O = (6*, 0) we have introduced the notation 57 = (6, —<fr) . 

Now we turn to the choice of our coordinate system in the overall center-of-momentum frame. We use the so-called 
"rigid body" system Kj s , associated with the final state plane spanned by the final three particles, in which the z-axis 
is taken to be the normal to this plane and parallel topx q\ . Thus the x- and y-axes are in the final scattering plane 
(see Fig. [I]). 

At a given three-particle invariant energy W, the relative orientation of the final particles within the final state 
plane is characterized by three independent variables for which we take the angle <f> p of the final nucleon momentum 
and the energies of the two mesons, U)\ and w 2 (see Fig. [T]). After straightforward algebra one obtains for the final 
nucleon momentum p 



p=\p\ = ^(W-cj 1 -u 2 y-M 2 N , (10) 
and for the relative momentum q of the two mesons 

q * = ±(u>f + u>l-M?-Ml)-^. (11) 

The orientation of the chosen coordinate system with respect to the beam axes may be specified by f2 7 = (</> 7 , 7 ), 
the spherical angles of the photon momentum k with respect to Kf s . One readily notes that in this coordinate system 
one has 8 p = 9 q = tt/2 and therefore 

D^ mp {<j> p -6 p ,-<j> p ) = (-iy- m rdl* mp (n/2) e -(-^, (12) 

D l ° m U q ,-e g ,-cf> g ) = (-!)-«£&_ (7r/2)e*"«*«. (13) 
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As will be shown soon, instead of <f> q only <j) qp = <f) q — <j) p is needed. It is related to uj\ and UJ2 by 

cos 4> qp = ^- (u% - ul - Ml + Ml ) , (14) 
2qp 

with p and q from Eqs. (ITU1) and (jTTJ) , respectively. Thus we will take as independent variables besides the photon 
angles f2 7 = (0 7 , </> 7 ) and <j) p the energies of the two mesons lo\ and UJ2 instead of p and <p qp and obtain the following 
representation of the T-matrix element making the angular dependence explicit 

T^^^u^He^-^e^ (15) 

JM 

with the contribution of the final partial wave 

= E ( J -v ) ( M - A A 4 ) ^-"^ 0^\l p3p m p ) , (16) 

which shows the explicit dependence on gp . Furthermore, we have introduced for convenience the notation 
O x A f J (l pJp m p ) = ( ~ 1) ^- JJ E i L (-l) l * +j * +l «%3 P %Ld l * mq (7r/2) 



X 



X ^^{{i^A)A\o^. (17) 



The following symmetry properties hold for the (l P jpm P ) 

M XLJ (l P J P m p ) = (-) L+l - +M - m "0^ J (l p j p m p ), (18) 
0^(l P j P -m p ) = (-y» +J O x A f J (l pJp m p ), (19) 

where the first one is a consequence of parity conservation. 

The symmetry relation of Eq. (|9]) leads to the following symmetry property of the amplitudes t 3 J£ 

t J -J-x-M qP ) = (-ir +M *^ M (-<M • (20) 

This means that for each J the number of independent amplitudes is 4(2 J + 1). 

The complex functions t^^f , depending on the meson energies u>\ and cl>2 only, provide a complete description of 
the process in a manner analogous to the description of a single meson photoproduction in terms of multipoles. It is 
worth to point out, that in contrast to the binary reactions the partial amplitudes are functions of the cm. energies 
of the final particles and, therefore, are to be determined for every point of the Dalitz plot. 

B. The differential cross section 

For the unpolarized differential cross section one obtains with the T-matrix of Eq. (|15[) 

dujidLU2dcosa~ ( ad>~ /v 4 Z — ' 

= /^jmi^lt^Yj^Oj^jp) (21) 
jm 

where we have defined 

s^ciw)^^) = ^c(w)3 £ i-iT M '( J M , _ J M J m ) 

J'M'JM ^ ' 

xEt-OWA-'/i fj, — X q)*^'^.^)**^!,^), (22) 
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with 



Ml 



C ^ = 4(2,)^ -Mj») ™ 

as a kinematical factor. One should note that the differential cross section depends on the relative angle (j> lp only 
besides on u)\, LJ2, and # 7 as is immediately evident in the absence of polarization effects. 
In terms of the electromagnetic multipole contributions one finds 



)<p q 

"p "OMpV" / "I " 



s i _ i y/ + J+M'+L'+L' 



X 

l 'pj' P m ' P l pjp m P 

E 

J'M'JML'L 

L V j 

A 





J 






L' 


1} 



E(~) A A -A m/ J '(O'X)* m- / (Wp^p)- (24) 



If with respect to the fixed final state plane only the direction of the final nucleon is detected, one obtains a semi- 
inclusive differential cross section by integrating the expression in Eq. ()21|) over loi and LU2 (setting without loss of 
generality (p p — 0, which means that 7 is measured relative to the direction of the nucleon momentum) 



d 4 a 



= S Sj m Yj m (^ 7 ) (25) 

jm 

as an expansion in terms of spherical harmonics in f2 7 with 



Sj m = / dwidw 2 S. jm (uJi, u 2 ) 

- ^c(w)3 E (- 1 ) 



-I M 

J'M'JM 



J' J j 

M' -M -m 



l p3p m p l p3p m P 
' + .T+M'+L'+L' 





J 




Hi 


L' 


1} 



x^(-l)^ r J ' m M J A j) / dwidwa^lf ("i^rt^i.a*), (26) 
or in terms of the multipoles 

^ = f ^ E E C-i/--*-* ■§•) (5 _X -tj {| j; f 

J p M p IL.iLTni.l'oj-ovn^ J 

x E (- 1 )' 7 " 

J'M'JML'L 

x E(-) A (a -X 0) / dwidw2 e i(m+M ^- 0^ ,J, (^X)* °m J (W^p) • (27) 

Since d 2 OQ/dQ. 1 is a real quantity, one has the property 

Sj m = {-) m Sj- m . (28) 

Furthermore, the cross section should be invariant under the simultaneous inversion of k and p, i.e. under the trans- 
formation 6 1P — > 7r — # 7p . Thus one finds as additional symmetry property 

Sjrti — ( )^ Sjm , (29) 
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from which the selection rule SV 



for j 



m 



- odd follows. This property can also be shown straightforwardly 



using Eq. 
namely 



(p?S|) with the help of Eq. ([!§]) . For identical mesons, one finds from Eq. (|2T))) an additional symmetry, 



Sj- T 



(30) 



which leads in conjunction with Eq. (|28j) to 3m Sj 



0. 



It is more convenient to use instead of the differential cross section the corresponding normalized quantity 

1 d 2 a 1 7 



W^y, - ... 

(To ally 

where the total cross section <tq is given by 



47T 



E 

jm,j>l,j-\-m—even 



(31) 



cr = 2v / 7rS'oo 

= Trc(^) ydwxdwa ]T _L_ |^(a, 1)W2 )| s 



(32) 



v\\i.JM 



and the expansion coefficients by 



77 



c(W) / dwidw 2 



vXfiJ'M'JM 



(33) 



Using the spherical harmonics expansion pip should enable one to interpret the experimental results without 
resorting to a particular model. This expression is an analogue to the expansion of the single meson photoproduction 
cross section in terms of Legendre polynomials. The coefficients Wj m are hermitesch functionals of the partial 
amplitudes i^*^. They obviously contain the whole information on the dynamics of the reaction with unpolarized 
particles and their values may in principle be extracted from the measurements and compared with model predictions. 
The selection rule Wj m = for j + m — odd may be used for a model independent partial wave analysis in the low 
energy region of the reaction, where usually only the first few waves contribute. 

Otherwise an integration over the angles # 7 and <pj gives the distribution of the events over the Dalitz plot 



d 2 a 



nc(W) 



T — 



jlt^i.wa)!' 



(34) 



Thus, as is well known, the partial waves of different J do not interfere in the Dalitz plot. In spite of its simplicity 
the expression in Eq. (|34l) can hardly be very useful in reconstructing even the modulae of the amplitudes Its 
use implies that one is able to establish a correspondence between variation of the amplitude as function of (wi,^) 
and a specific value of the total angular momentum J. Obviously, for this purpose a detailed model is needed which 
relates J to particular decay channels. In this sense, using the moments Wj m should be more promising. 

It is also clear that the information on the unpolarized differential cross section only is insufficient for a model 
independent determination of the amplitudes tjjjjf. In the general case of photoproduction of two pseudoscalars eight 
independent complex functions are required to fix the spin structure of the amplitudes. Since the overall phase is 
always arbitrary, one has to measure 15 independent observables at each kinematical point. However, in certain cases, 
e.g., when the reaction is dominated by a single partial wave, using the moments Wj m enables one at least to draw a 
qualitative conclusion with respect to the partial wave structure. As an illustration, we consider in the next section 
the theoretically interesting case of 7r 7r° and ir°r] photoproduction on a proton. 



III. APPLICATION TO 7 p ->■ TrVp AND 7 p -> ty°t]p 

The measured total cross section for jp —> 7r°7r°p exhibits a rather steep rise in the energy region below the 
Di3(1520) resonance (see, e.g., 23]). At the same time, the existing models with a dominant contribution from 
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FIG. 2: The moments Wj m for "/p — > n ir p as functions of the photon lab energy, normalized such that Woo = 1. 



Di3(1520) and a moderate role of the Roper resonance predict a cross section which increases rather slowly with 
increasing energy and is, therefore, far below the data. It is reasonable to assume that the almost linear energy 
dependence of the data indicates a contribution of a large fraction of s waves in the final state. The main mechanism 
providing the s-wave part in tttt photoproduction is the A Kroll-Ruderman term, appearing after minimal substitution 
of the electromagnetic interaction into the itNA vertex. This term, however, vanishes in the neutral channel. The 
situation is similar to that in single ir° photoproduction at low energies. Here the Kroll-Ruderman does not enter the 
amplitude, thus leading to a visible suppression of the cross section for 7p — > ir°p in comparison to the ir + or ir~ case. 

A possible large contribution of the Roper resonance Pn(1440) in the region E 1 = 500 — 600 MeV as assumed in 
Ref. [ll| seems to be excluded by more resent analyses. Furthermore, this assumption should be in disagreement with 
the experimental results of Ref. [24| for the helicity dependent total cross section Act = ct 3 / 2 — cri/ 2 - There it was 
found that in the energy region up to at least E~ ( — 800 MeV the 3/2 part dominates over the 1/2 part. This means, 
that the P\\ wave, which contributes only to (J 1 / 2 , should be overwhelmed by the waves with higher spins. 

Thus the question concerning the partial wave structure of the amplitude for — > ir°ir°N is still open. In order to 
reveal in this case the mechanism responsible for an unusually large fraction of the s wave part in the it it amplitude, 
it is useful to analyse the moments Wj m throughout the energy range from threshold up to the -Di 3 (1520) peak. In 
order to keep the number of parameters limited, one can use only the lowest partial waves. Their choice is inspired 
by the previous isobar model analyses of Refs. [13, [H, [l4| showing that only waves with J < 5/2 are important below 
£ 7 = 1 GeV. 

As an example we show in Fig. [2] the variation of Wj m for j < 3 as predicted by the tttt model of Ref. [14|. The 
model [IH is based on a traditional phenomenological Lagrangean approach with Born and resonance amplitudes 
calculated on the tree level. The interaction within the nN and 7T7t pairs is effectively taken into account via A, p and 
ct. The ttttA^ state is then produced through intermediate formation of 7rA, pN and ctA^ channels. The contributions 
from the resonances are parametrized in the usual way in terms of a Breit-Wigner ansatz with energy-dependent 
widths. For the parameters of the model, i.e. masses, partial widths and electromagnetic couplings of resonances, the 
corresponding average values from the compilation of the Particle Data Group were used. 

In case of it it production due to the identity of the two mesons we have an additional symmetry relation 



W{6- I ,<t>- I ) = W{6^2TT-4> 1 ), 



(35) 
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FIG. 3: Same as in Fig.[2]for 7p — ¥ TT°T]p. The dashed lines represent the imaginary parts. 



which is a consequence of the symmetry property in Eq. (|30[) . The moments for j = 3, 4 are small as are those 
for higher values of j which are not shown. In the region E = 650 — 800 MeV the moments Wn and W20 exhibit a 
crucial energy dependence due to the -Di3(1520) resonance, dominating the reaction 7p — > 7r°7r°p at this energy. Large 
values of the moments with j odd indicate the presence of waves with opposite parities. In particular, the structure 
in Wn is due to an interference between the wave J p = 3/2~ dominated by Di3(1520) and the waves J p = 1/2+ 
and 3/2 + . The latter are saturated, apart from the Roper resonance, by the Born terms. The contribution of Wn 
becomes minimal in magnitude in the region around E 1 = 650 MeV, where the real part of the -Di3(1520) propagator 
vanishes, and it interferes weakly with the predominantly real Born amplitudes. Thus, if our notion about the 7r°7r° 
photoproduction mechanism is correct we expect a rather small value of the moments W20 and W22 and a relatively 
large value of W\\ in the region below the Di3(1520) peak. 

In this respect we would like to note that according to the fit in Ref. [16| there must be a large contribution of the 
resonance £>33(1700) to the channel tt°tt°p in a wide energy range from the lowest energies up to E y — 1.4 GeV. In 
particular, inclusion of this resonance into the amplitude explains both the steep rise of the total cross section below 
.Ely = 700 MeV and the second peak observed at E~ = 1.1 GeV. If the resonance 1)33(1700) is indeed so important in 
the 7r°7r° channel, it should increase the values of W20 and W22- All in all, a measurement of these moments will help 
us to understand the role of d-wave resonances with J = 3/2 in 7r°7r° photoproduction. 

As for tt°v photoproduction, the partial wave structure of the corresponding amplitude was investigated in detail 
in Refs. @, EEEl- There it was shown that the J p = 3/2~ wave, containing D 33 (1700) and probably £> 3 3(1940), 
apparently dominates the reaction in a wide region from threshold to about E 7 = 1.7 GeV. Other waves, primarily 
l/2 + and 5/2 + , manifest themselves in angular distributions of the final particles mostly via interference with the 
dominant 3/2" wave. 

In Fig. [3] we present the energy dependence of the expansion coefficients for jp — > ir°r]p obtained using the isobar 



model of Ref. [17]. Here the relation (f3"0"f does not hold, so that the moments Wj m with m ^ have nonvanishing 
imaginary parts (dashed lines in Fig. [3]). The calculation follows the same line as for the 7r°7r° case. Namely, the final 
ir i]N state results from the two step decay of baryon resonances via the intermediate quasi-two-body channels rjA 
and 7r°S'ii(1535). The parameters of the model were fitted to the angular distributions of the final particles measured 
in Ref. [6j. The fitting procedure is described in [17| and the reader is referred to this work for more details. 

Firstly, as one can see in Fig. [3] in spite of the mentioned dominance of the 3/2~ wave, the values of W20 and W22 
are small. This is because of the closeness of the 3/2 and 1/2 helicity couplings of the resonance Di3(1700) (see, 
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e.g., the discussion in Ref. As a result, the hermitian forms of t%j^f entering W 2 o and W22 according to (f33|) 

almost cancel each other. At the same time, we obtain a rather large value of the coefficient Wu, mainly determined 
by the interference between the resonances D 33 (1700) and P 3 i(1750). According to these results we may expect that 
the data for ir°r] will show relatively small values of all moments except for W\\. If this prediction is not confirmed 
by measurements one has to critically review the existing conceptions about the dynamics of photoproduction, 
based on the results from Refs. @, 0, E E3, d • 

IV. CONCLUSION 

Practical methods for the analysis of the partial wave structure of reactions with three particles in the final state are 
obviously needed for the study of the dynamical features of two-meson photoproduction. The formalism used in the 
present paper specifies the final ttttN states by means of two cm. energies and two angles, determining the orientation 
of the final state momentum triangle (final state three-particle plane) with respect to the beam axis. The partial wave 
decomposition may then be performed via a transition from the continuum variables (angles) to the set of discrete 
variables JM being the total angular momentum J and its projection M on the normal to the three-particle plane. 
The corresponding partial wave amplitudes i^*^ contain the whole information on the production dynamics. We 
would like to stress the fact that this method does not involve a decomposition with respect to the angular momenta 
of the final two-body subsystems and is in principle free from any assumptions about the production mechanism. 

In the present paper we have considered only the unpolarized differential cross section. Although this quantity 
does not allow a unique determination of the amplitudes t'l^, the information on the angular distribution of the 
participating particles can serve to place restrictions on contributions of states with definite angular momentum and 
parity. This in turn is crucial for our understanding of the resonance content of the reaction. For this purpose the 
differential cross section has been expanded in terms of spherical harmonics with coefficients or moments Wj m in a 
manner similar to the representation of the binary cross section in terms of Legendre polynomials. 
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